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We argue that the Milky Way (MW) contains thousands intermediate mass black holes (IMBHs) and mini- 
halos with a fraction of IMBHs still being enshrouded in extremely dense mini-spikes of dark matter (DM) 
particles. Each containing lO^M© of dark particles and no baryons in a sphere of 50 pc radius, the minihalos 
are dense enough to survive the Milky Way tide with the nearest minihalo within 2 kpc from the Sun. The 
IMBH is formed off-centre in a minihalo by gas accretion, and its growth adiabatically compresses a finite 
density of surrounding dark matter into a r^^'^ mini-spike. Some IMBHs recentre on their minihalos after 
dynamical friction, and some IMBHs are ejected by birth kicks. Detectable by GLAST, the mini-spikes and 
minihalos should stand out the background and dominate the neutralino annihilation in the smooth MW and 
satellite galaxies. If they are the unidentified EGRET sources, upper limits can be set on the branching ratio of 
neutralino annihilation. The supermassive BH of the MW, if orginates from an IMBH, is also likely enshrouded 
with a mini-spike. 

PACS numbers: 95.35.+d: 98.35.Gi; 98.35.Jk 



Supermassive black holes (SMBH) likely grow from 
smaller intermediate mass black holes (IMBH) by gas ac- 
cretion; to be massive enough and early enough, a seed of 
^ lOOOM© at z ~ 20 is necessary to grow to a 3 x lO^M^ 
SMBH at z = 6 via Eddington accretion of gas (e.g., Qll^ 
and references therein). Gas accretion is favoured because 
of the empirical result [^j^ that requires radiatively efficient 
accretion in order to account for both the quasar luminosity 
function and the mass in black holes at present. Conversely, 
once the quasar/ AGN phase is over, black holes undergo ra- 
diatively inefficient accretion as is the case for Sgr A* (e.g., 
fs*]). Here we explore implications of growth of SMBH by 
accretion in the protogalactic environment. We will study the 
evolution of the seed IMBHs. These smaller IMBHs are plau- 
sibly formed in a subset of the first objects, where accretion 
rates are high. Hierarchical merging of minihalos leads to a 
population of IMBHs in MW-sized halos. A mini-spike of 
CDM is inevitable around the IMBH due to adiabatic growth 
in the minihalo. The mini-spike is robust to halo merging, and 
survives around IMBHs that populate the halo. The subse- 
quent role of BH merging depends on dynamical friction and 
loss cone refilling. These processes are likely to nurture the 
growth of the central seed IMBH into a supermassive BH. 

The IMBHs, and in particular, the dark matter spikes, are 
potentially observable via neutralino annihilations. We men- 
tion two specific examples. Firstly, a central IMBH is inferred 
to be associated with a group of massive stars IRS 13 offset 
by about 2 arcsec from the Galactic centre (SgrA*) 0|. Sec- 
ondly, the Whipple/HESS detection of > 3TeV gamma rays 
from the Galactic centre (|7, 8, 9]) would require a high an- 
nihilation cross section ^^10 pb if this emission is due to a 
neutralino annihilation signal produced by an halo 
profile 1 1 1 ] . Unfortunately, to create the ACDM relic density 
requires a cross-section of < Ipb. Nor even for subdomi- 
nant neutralinos can SUSY theory easily produce such a large 
cross-section above 10 TeV. A spike may provide the best way 
of reconciling the high annihilation flux with the dark matter 
relic abundance lflllfl2ll . 



The redshift and mass of the first clouds, stars and IMBHs: 
The temperature T of the first clouds that contain gravitation- 
ally bound baryons depends mostly on the virial mass M with 

AI = 1O^M0 , where we have scaled with typ- 

ical values of the redshift z and temperature T; cooling by 
H2 requires a virial temperature in excess of several hun- 
dred degrees K. The collapse redshift for clouds of mass AI is 
ss O.Giy X 16(A//10^Mq)^ where e < 1 and i/ is the height 
of the gaussian peak. Hence the desired redshift is sensitive 
primarily to the peak rarity, and must be in the range 10-30. 
The WMAP-inf erred optical depth of the universe to Thom- 
son scattering also suggests a redshift of 15-20 for the first 
objects that reionised the universe. 

The mass of the first objects can be estimated by noting that 
the accretion rate is of order v^/G ^ 10"'* — lO^'^Mo/yr, 
where we take Vs to be the virial velocity. Accretion can only 
be halted by feedback, and it takes at least 10® years to gen- 
erate feedback from massive stars, hence the masses of the 
first stars are of order 100 — IOOOMq, as found in numeri- 
cal simulations of primordial gas clouds by 1 13]. Finally, in 
clouds of lO'^K where Lyman alpha cooling dominates, the 
accretion rate onto the protostellar core is around O.OSMo/yr, 
hence feedback cannot inhibit the formation of an IMBH, of 
mass up to ^ 1O''M0. 

Annihilation flux of (mini-}halo models with finite core and 
plateau: The (bolometric) annihilation flux integrated over a 
sphere of WIMPs of mass and number density n{r) = 
at radius r is 



L^^{< r) — 'inr'^dr (crv) m^c^ 
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where annihilation imposes a plateau of size r < on the 
present-day dark matter density because neutralino number 
density above = ^ '""^^S"^'' 10«cm"3 will 
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annihilate away over a Hubble time (e.g., I14i llSll ). where 
the cross-section (crv) /c Ipb, fixed by the relic abundance 
fi^ required in a ACDM cosmology. The total detectable flux 
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quickly converges at large radius since p{r) generally drops 
as a steeper than r 
typical SUSY mass 



^■^ power-law beyond Vp. Adopting the 
= 50GeV, we have 
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Note that less than 10% of the above calculated flux for neu- 
tralino models will likely go into photons (7-rays). More gen- 
erally and crudely speaking, L^^{< r) oc 

We model the density p{r) of the minihalos and the MW 
halo with a subset of models of Zhao|16] (also known as 
generalized-NFW profiles). To be conservative in estimat- 
ing the annihilation signal, we choose a subset of initial halo 
models with a finite core and without any outer truncation. 
Our models are also inspired by the recent exponential-like 
model of Navarro et al. (| 17], N04 model hereafter), but are 
easier to work with analytically. Specifically, we model halos 
with density at radius r given by 



p{r) 



[x/2 + 1/2] 
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Here Vc is, as in an NFW model, where the local density 
power-law index d log p/d log r — ~2. Our models are speci- 
fied by two parameters Vc and pc = p{rc), and the virial radius 
Tvir = cVc is given implicitly by 



3Mvir 



176pcrit(2) 



2p(rc)(l + c3/5 + c3/55) 



i.e., the mean density p(r 
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is approximately 176 
times the critical density PcritC-z) of the universe at a given for- 
mation redshiftz. Note the density ratio p(rc)//o(rc) = 2.436, 
2.3178 and 2.4644 for our cored model, NFW model| 10], and 
N04 model. Unlike the NFW model, our model has an outer 
density profile without truncation, and a finite central den- 
sity and a finite asymptotic mass. For models with the same 
core radius and density nomalisation p{rc), A/(rvii) differs 
by only 10 percent between NFW models and our model for 
typical concentration c = 1 — 1000. 

There are a few advantages of our model. It is straighfor- 
ward to derive and show that the potential, the circular veloc- 
ity curve and even the velocity dispersion curve of this class 
of halo models are all simple analytical rational functions of 
r^/'^, while some special functions enter the NFW and N04 
models. As shown in Fig.l, our model reproduces the density, 
mass, velocities of NFW and N04 very well over the range 
O.Olrc < 7' < lOOrc at least for systems with c > 1. 

Formation of minihalos and massive black holes: Current 
numerical simulations suggest that the first stars form from 
H2 clouds in dark minihalos of mass about 3 x 10^/i~^Mq, 
virial radius of 60ft-^^pc and virial velocity of 5 kms^^, 
and initial virial temperature of 1000 Kelvin (for a molecu- 
lar weight of unity) at a redshift of z ~ 20, when the uni- 
verse was about 0.15 Gyr old. The cooling time is compa- 
rable to the dynamical time, both being of order 1 Myr, and 
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FIG. 1: log-log radial profiles of the rescaled halo density 
p{r) / p{rc), escape velocity Vcscir) / {K\/2), dispersion 2(j{r)/K, 
circular velocity Vciy:c{r) / K , and enclosed mass M{r) /lQQp(rc)r't 
in our cored halo model (crosses), NFW model (blue dashed lines) 
and N04 model (red lines). Here K = Gp(rc)rc, and Tc is the 
core radius. 



self-gravitating isothermal clouds of order a thousand solar 
masses are able to develop. Once the H2 metal-free gas in 
these systems cools to a temperature of ~ 200K and density 
of ^ lO'^cm^'^, the H2 cooling saturates and subfragmenta- 
tion ceases L18il . These clouds are generally not at the exact 
centres of the minihalos, but have dynamically decayed to the 
central parsec of the minihalos with typical random motions 
of less than 10 kms^^. The first stars are generally massive 
(about 100 — IOOOM0) because of inefficient cooling with 
zero metals. They are short-lived. According to ] 19] (see their 
Fig. 2), these very massive zero-metalicity stars have large he- 
lium cores, which are unstable to pulsations due to electron- 
position pair instability. A supernovae leaving no remnant 
happens for 140 — 260Mq stars, and a more massive star col- 
lapses into an IMBH after lossing about half of the initial mass 
by wind. However, Lyman alpha-cooling clouds would be ex- 
pected to form more massive IMBHs, up to 10'^ — 1O^M0. 

There are several plausible outcomes for the IMBH and 
the gas inside a minihalo. In some cases, the IMBH could 
be ejected out of the shallow potential well of the minihalo 
if born with a significant kick/recoil velocity (greater than 
the 30 kms^^ escape velocity at the centre of the mini- 
halo) L2Qil . In other cases, the IMBH spirals to the cen- 
tre of the mini-halo due to dynamical friction (the timescale 
tdf ~ (t/pc)'^ ~ IMyr for a thousand solar mass BH formed 
at Ipc). At the centre or en route, these IMBHs may accrete 
gas at a super-Eddington rate with a mass-doubling time scale 
400e ln(2) sa 1 — 20Myr for a plausible range of the accretion 
efficiency e — 0.005 — 0.1. Eventually after some uncertain 
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period of order 10 Myr, the first SNe progenitors are born in 
these minihalos, and shortly afterwards (about 4 Myr) most 
of the gas in the minihalos (perhaps with the exception of the 
tightly bound accretion disk around the IMBH) is cleared out 
when the SNe explode. Thus the only remnants of the first 
generation of star formation are a population of minihalos and 
IMBHs, which are either well-centred due to dynamical fric- 
tion or well-separated due to strong kicks. To be specific, 
we consider minihalos specified by eq. (|3j with rc — 1.6pc 
and p{rc) = 500Mopc^^ (comparable to an NFW halo of 
Tvir = 50 X Tc = 80pc). We assume conservatively that 
10 percent of the minihalos harbour IMBHs which eventually 
grow to a final mass of IQ-^Mq or IO'^Mq. These values are 
plausible for a rm.s. kick velocity ~ 100 kms~^, and accre- 
tion efficiency e ^ 0.02 in the mini-halo in 10 Myrs. So there 
are of order 10'^ minihalos containing IMBHs. 

Adiabatic compression of dark matter in minihalos: Now 
consider the effect on the dark matter. The formation and the 
growth of the IMBH will compress the surrounding dark mat- 
ter. The compression is close to adiabatic since the forma- 
tion of the IMBH and the subsequent accretion are on time 
scales of the cooling time and Salpeter time; both are some- 
what longer than the dynamical time. Given the finite density 
and phase space density of dark matter at an off-centre posi- 
tion, we expect a moderate r^'^/^ power-law dark matter den- 
sity distribution surrounding the IMBH. This type of cusp (or 
"spike") is the result of growth in a uniform background! 25 1; 
a steeper cusp occurs if a well-centred BH grows adiabati- 
cally (e.g., |151). Such a mini-spike will survive dynamical 
friction, and be added to the central core of the minihalo once 
the IMBH recenters. We expect of order of m, ~ IOOOM0 
of dark matter in this spike. So the dark matter mass profile 
with a mini-spike is 

M'^{t) =(^l + M(r), M{h) = m., (5) 

where b = 0.045pc or 0.14pc so that M"p(6) 2M{h) = 
m. = IO^Mq or IO^Mq. 

We consider several possibilities for the dark matter Fig. 2 
shows the resulting dark halo density and the annihilation flux. 
In the simplest case, the minihalos have the pristine cored den- 
sity profile; the IMBH either never formed or has been kicked 
out of the system. For a minihalo which keeps its IMBH, the 
IMBH draws in a mini-spike of dark matter with a plateau 
inside about 0.005pc depending on its mass. The flux per log 
radius bin is of order lO'^L©, and most of it comes from within 
O.Olpc of the minihalos. 

Observability of minihalos inside the MW: By a redshift 
of z — 1 — 6, nearly all the minihalos and IMBHs are 
bound inside galaxy-sized objects. About 10'^ — 10^ should 
be bound inside a Milky- Way-sized galaxy |21], which we as- 
sume to have a density given by eq. (|3} with rc — 11 kpc and 
p(rc) = O.OO55M0pc^^ (this is consistent with the local DM 
density and resembles an NFW model of r^ir = 8rc). These 
minihalos should be distributed in a way similar to the overall 
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FIG. 2: Left: log-log plot of the radial profiles of spherically 
averaged density (black lines) in Mopc"^ inside a mini-halo in 
three cases: a strong mini-spike around a lOOOM© IMBH (thick 
line), or a moderate mini-spike around a IOOM0 IMBH, or a pure 
cored minihalo where the IMBH either was never formed or was 
ejected (thin line). Lines (green parallel) of constant enclosed mass 
IM0, IOM0, W^^Mq are shown to guide the eye. For compari- 
son are the average densities inside a dwarf galaxy (lower thick line, 
with dispersion of 10 kms~^ and a core of 200pc) and inside the 
pristine MW halo (lower right thin line). Right: the logarithm of 
the reseated annihilation (bolometric) luminosity x gjg^y 

(nearly horizontal red lines to the right) of these models. Also shown 
are the total luminosity of 10^ minihalos with IOOMq IMBHs (thick- 
est black at top middle), the luminosity of the MW as a whole (thick 
red) or scaled down by a factor for a cone of f2 = 10~^ (thin 
red). The minihalos stand out as clumps above the MW background. 



dark matter profile of the galaxy, producing a clumpy anni- 
hilation brightness distribution; we assume the stars and the 
central SMBH in the MW form later. 

A typical nearby minihalo at about 3 kpc from the Sun 
sustains a solid angle Ail ^ 10~^ (about the resolution of 
EGRET) comparable to the angular size of dwarf galaxy satel- 
lites at lOOkpc, but minihalos are much more luminous and 
brighter in annihilation than dwarf galaxies because of mini- 
halo's much denser core than dwarf galaxies; the latter have 
central dark matter densities of (0.05 — 0.5)Mqpc^'^ |26]. Ex- 
cept towards the Galactic center a minihalo stands out in den- 
sity and annihilation flux from the MW background in a beam 
of Aft < 10^'^ at all distances (cf. Fig. 2). A population of 
10^ minihalos with lOOM© IMBHs is about 100 times more 
luminous than the MW background as a whole. These mini- 
halos are low-mass (^ IO^Mq) but dense (about IMopc^"^ 
on a scale of 30pc) enough so that 99 percent of the minihalos 
survive intact in the halo, and follow the overall distribution 
of dark matter. There are of order 10 minihalos within 1 kpc 
of the MW centre, where dynamical friction could bring one 
of these minihalos to the centre of the MW within a Hubble 
time (cf. Fig. 3 of L15J ). without being slowed down by tidal 
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stripping '221; note that these minihalos are typically denser 
than the MW halo on a scale of 60pc. If the IMBH of the 
minihalo could serve as a seed for the supermassive BH in the 
MW. 

The isolated IMBHs are virtually invisible. The X- 
ray luminosity due to Bondi-Hoyle acretion of an IMBH 
passing through the MW gaseous disk is estimated to be 
(e/O.l)(n/3Ocm-3)(r7i./3OOM0)2L0 |22il. So the signal is 
much fainter than the million-solar-mass BHs of Lacey & 
Ostriker|231. Such a low luminosity is just detectable with 
present facilities if the IMBH is close-by (say 2 kpc), but the 
chance of a BH passing through a cold H2 or HI cloud (with 
filling factors 0.3% and 2% respectively) at a given time is 
very low. The local density of IMBHs is about 0.3 kpc^'^ as- 
suming 10* IMBHs in the whole halo, and we expect only of 
order one IMBH hitting a cloud with a filling factor of '--^ 1 % 
in the entire disk of 8 kpc in radius and 200pc in thickness, an 
approximation for the cold H2 and HI cloud distribution. Also 
being much fewer than the million-solar-mass BHs of li23ll . the 
minihalos do no damage to the MW disk. 

Despite being a dynamically subdominant population, 
the minihalos and IMBHs speed up two-body relax- 
ation in the MW. The relaxation time is given bv 12311 

IMyr _ P(r)<y-'' M_ ^Up,p A> _ 

trlx ~ lO6M0pc-3(lOOkms-i)-3 IO^Mq' "'"'^re IVl — 

M^+r.' ~ 1 - IOOOMq is the mass density-weighted 
mean mass of deflectors depending on tidal stripping; the 
lower value is for minihalos well inside 1-10 kpc of the Galac- 
tic centre, and Jm ^ 10^'^ is the fraction of matter mass 
density in the form of minihalos of mass M . The heav- 
ier the deflectors, the faster the relaxation. So compared to 
the case of purely stellar deflectors, the relaxation is about 
M/M0 ~ 1 — 1000 times more efficient. This might help 
to refill the loss-cone of the central BH by scattering stars and 
dark matter particles! 24 1 . To refill the loss-cone of 3 x lO^Mo 
BH on the scale of the Ipc would require a relaxation time of 
order a Hubble time. 

Our model predicts around 1500 L0 (bolometric) per so- 
lar mass in spike neutralinos. Observationally, the annihila- 
tion signal is best seen via gamma rays. EGRET observa- 
tions set upper limits on unidentified sources, in particular on 
the source near the Galactic Centre, in gamma rays above 1 
Gev of around IOOLq iIztIi . For our models to be consistent 
with observational data in the MW and a gamma ray branch- 
ing ratio of order 10 percent, we require of order a solar mass 
of the neutralinos in the spike for typical SUSY neutralinos 



with m. 



50GeV. This is consistent with our models of 



adiabatic growth of IO^Mq IMBHs and models with ejected 
IMBH. The data would also be consistent with the IO'^Mq 



IMBH model in case of heavier neutralinos (m. 



lOTev 



| 28|). Evidently the GLAST satellite, with two orders of mag- 
nitude more sensitivity than EGRET, will be capable of setting 
stringent constraints on halo IMBH. 
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